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Two-dimensional (2D) crystalline colloidal monolayers sliding over a laser-induced optical lat¬ 
tice providing the periodic “corrugation” potential recently emerged as a new tool for the study of 
friction between ideal crystal surfaces. Here we focus in particular on static friction, the minimal 
sliding force necessary to depin one lattice from the other. If the colloid and the optical lattices 
are mutually commensurate, the colloid sliding is always pinned by static friction; but when they 
are incommensurate the presence or absence of pinning can be expected to depend upon the system 
parameters, like in one-dimensional (ID) systems. If a 2D analogy to the mathematically estab¬ 
lished Aubry transition of one-dimensional systems were to hold, an increasing periodic corrugation 
strength Uq should turn an initially free-sliding, superlubric colloid into a pinned state, where the 
static friction force goes from zero to finite through a well-defined dynamical phase transition. We 
address this problem by the simulated sliding of a realistic model 2D colloidal lattice, confirming 
the existence of a clear and sharp snperlubric-pinned transition for increasing corrugation strength. 

Unlike the ID Aubry transition which is continnous, the 2D transition exhibits a definite first-order 
character, with a jump of static friction. With no change of symmetry, the transition entails a 
structural character, with a sudden increase of the colloid-colloid interaction energy, accompanied 
by a compensating downward jump of the colloid-corrugation energy. The transition value for the 
corrugation amplitude Uq depends upon the misalignment angle 9 between the optical and the 
colloidal lattices, superlubricity surviving until larger corrugations for angles away from the ener¬ 
getically favored orientation, which is itself generally slightly misaligned, as shown in recent work. 

The observability of the superlubric-pinned colloid transition is proposed and discussed. 
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I. INTRODUCTION 

The great progress of nanofriction of the last decades 
has enabled increasing insight into the microscopic nature 
of frictioniir— In particular, the sliding between contact¬ 
ing crystalline surfaces has become of higher interest af¬ 
ter the manipulation of nanoscale-sized objects revealed 
how the atomistic features of an interface may play the 
ultimate role in determining its tribological properties. 
Fresh interest is now generated by artificial but highly 
instructive systems represented by 2D colloidal monolay¬ 
ers sliding in an optical lattice, the main system which we 
use as a working example in this study4i Two atomically 
flat and perfect crystal surfaces in contact are commen¬ 
surate when the mismatch ratio p between their lattice 
parameters is rational. Perfect matching (p = 1) between 
two lattices of same symmetry represents the simplest sit¬ 
uation, where the atoms of one crystal lattice (the slider) 
perfectly fit the minima of the two-dimensional “corru¬ 
gation” periodic potential landscape generated by the 
other lattice (see Fig. [T^). In this case a finite and large 
external force is required to dislodge the slider atoms 
from the potential minima, so as to nucleate, at finite 
temperature)^!^ the onset of sliding. 

A different situation arises when the periodicities of 
the two contacting lattices are mismatched, as is realized 
either by a value p 1, or else by a mutual rotation 
of two lattices, as shown schematically in Fig. [Ud. The 
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FIG. 1. (Color online) (a) Schematics of a p = 1 commensu¬ 
rate interface between a colloidal monolayer (black dots) and 
a triangular optical lattice potential (light spots correspond 
to the minima), (b) A colloidal monolayer of lattice spac¬ 
ing Oc with a mismatched triangular potential of periodicity 
ai = puc. The two lattices are mismatched but aligned, (c) 
Colloidal monolayer with a small misfit angle 9 with respect 
to the mismatched substrate. This kind of rotation generally 
lowers the energy and is a likely equilibrium feature.— 


present study concerns interfaces that are incommensu¬ 
rate. A situation far more common in nature, it also 
naturally exhibits more interesting tribological features 
than commensurate ones?^ We focus especially on static 
friction, the minimal sliding force necessary to depin one 
lattice from the other. The main issue we address here is 
the nature of static friction against sliding of a 2D lattice 
onto a 2D incommensurate corrugation. 

The reference point for our understanding is the ID 
case, epitomized by the so-called Frenkel-Kontorova (FK) 
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model, consisting of a harmonic chain of classical point 
particles in a static sinusoidal potential, an idealized sys¬ 
tem whose physics has been very widely studied and 
understood)^ In ID, the incommensurability between in¬ 
terparticle spacing and sinusoidal potential wavelength 
gives rise to a potential-induced deformation of the chain 
particle position, that can be described by a deformation 
of the chain’s phase $(a:) relative to that of the poten¬ 
tial, from the non-interacting straight line of slope {p—1) 
to a staircase with the same mean slope but formed by 
nearly horizontal, approximately constant phase terraces. 
These commensurate domains are separated by steps in 
the relative phase, often called solitons (or kinks) and an- 
tisolitons (or antikinks), where most of the misfit stress 
is concentrated. 

This ID system is known to display a remarkable 
dynamical phase transition, first described by Aubryi^ 
When the corrugation is weak, below a critical magni¬ 
tude, when the soliton widths are comparable to or larger 
than that of terraces and the overall phase modulation 
is small, the FK chain exhibits no pinning, the static 
friction is zero. In this regime one finds particles at all 
potential values, and all possible phases of the chain par¬ 
ticles are accessible in dynamics. Since the total poten¬ 
tial energy is independent of the phase, 5E/5^ = 0, the 
chain as a whole can be displaced by any infinitesimal 
force. This absence of static friction has been dubbed 
“superlubricity” However, above a critical corrugation 
magnitude, whose value depends on precise parameters 
and incommensurability, the chain develops a nonzero 
static friction through a continuous “Aubry” phase tran¬ 
sition. The essence of the Aubry transition is that the 
probability to find a particle exactly at a potential maxi¬ 
mum drops mathematically to zero (of course at T = 0). 
Thus, even if 5E is still zero, this new dynamical con¬ 
straint limits phase-space accessibility, breaking ergodic- 
ity in a way similar to a glass transition, with the onset of 
static friction and of chain pinning against free sliding. 
This transition is of a dynamical nature, so that even 
if the step-terrace deformation of the chain’s phase <I>(a;) 
becomes more marked in the pinned state, it still remains 
similar in nature to the superlubric one. Although sug¬ 
gestive, the static structure difference between the two 
does not reveal the essence of the transition. A more rel¬ 
evant signature is provided by the sliding friction. In the 
superlubric state, sliding of the chain occurs under any 
applied force however small, and can be envisaged as a 
state of flow of the solitons. Soliton motion causes dissi¬ 
pation, increasing with speed in correspondence with in¬ 
creasing emission of phonons. In the pinned state, chain 
sliding occurs only once the static friction force is over¬ 
come, at which point the sliding dynamics becomes gen¬ 
erally quite different from the superlubric state. 

In view of this well-known ID case, it is physically 
clear that a qualitatively similar state of affairs should 
also occur in 2D incommensurate sliders. That is actu¬ 
ally a case much more relevant in practice, where one ex¬ 
pects free, superlubric sliding for weak corrugation, and 


pinning with stick-slip for strong corrugation. Superlu¬ 
bricity is documented in several real 2D systems. Very 
small values of the static friction force Eg were observed 
for a rotationally misaligned graphene flake sliding over 
graphite^ii the pinned configuration corresponding to the 
aligned and therefore commensurate case. Superlubricity 
and consequent ultra-low dynamic friction has also been 
demonstrated or implied in a number of cases, such as 
also telescopic sliding among carbon nanotubeS j^^d^ in 
cluster nanomanipulation studiesj^d^ and in the sliding 
of rare-gas adsorbatesd^iil However, there is no experi¬ 
mental case so far exhibiting a clear superlubric-pinned 
transition among 2D incommensurate sliders under fixed 
geometrical conditions. An experimental system of fresh 
relevance where that goal might be pursued is represented 
precisely by two-dimensional crystalline colloidal mono- 
layers sliding over a laser-induced periodic corrugation 
potential. Measurements^ and computer simulations^ of 
colloidal monolayers interacting with a laser-induced pe¬ 
riodic potential demonstrated the possibility to attain 
negligible values of the static friction force in incommen¬ 
surate geometries. Therefore a study of the transition 
from the superlubric to the pinned state, done at con¬ 
stant lattice mismatch, and at constant mutual align¬ 
ment is potentially feasible for sliding colloid monolayers. 
These systems are excellent candidates as they allow to 
change all the relevant parameters of the interface: the 
corrugation strength Uq and the colloidal crystal lattice 
parameter Oc, which sets the mismatch ratio p = oi/oc 
given the optical potential periodicity ai, fixed by the 
laser. The mutual alignment or misalignment between 
the colloid and the optical lattices is decided by energet¬ 
ics. As was established by a recent study, a small mis¬ 
alignment generally prevails, directly detectable by the 
colloid moire pattern?!^ This kind of misalignment, de¬ 
picted in Fig. [1]:, is well established in surface-adsorbed 
rare-gas monolayers 

The question we address here, only marginally touched 
upon so far— is whether at constant mismatch the weak- 
and the strong-corrugation regimes will again or will not 
be separated by a well-defined superlubric-pinned phase 
transition of the same type as that described by Aubry 
in ID - the only case where a mathematical treatment 
has been worked out. Colloid monolayers offers the ideal 
chance to verify this point if not analytically at least 
numerically in a realistic and relevant case. 

For a 2D triangular crystal monolayer interacting with 
a mismatched rigid corrugation periodic potential the ID 
soliton-terrace staircase is replaced by a 2D superlattice 
of nearly commensurate domains. The domains are sepa¬ 
rated by soliton (antisoliton) misfit dislocation lines, de¬ 
pending on the specihc value of p, and of the misalign¬ 
ment angle. In triangular symmetry, three families of 
soliton lines, oriented at ±120 degrees from one another, 
will generally coexist. The coincidence plot between the 
two lattices (monolayer and periodic corrugation) real¬ 
izes a moire pattern, consisting of a patchwork of nearly 
commensurate, roughly hexagonal domains separated by 


3 


the soliton lines, which concentrate the mismatch of the 
two lattices)^ Like in ID, the sliding of the 2D monolayer 
enacts a ’’flow” of the solitons in the moire pattern. A 
spontaneous angular misalignment, which when present 
lowers the energy by increasing the interdigitation of the 
two lattices, strongly modifies that pattern and corre¬ 
spondingly increases frictionJ^ 

The main questions which we plan to address here are: 
(a) Is there a sharp superlubric-pinned transition in the 
sliding of incommensurate colloid monolayers for increas¬ 
ing corrugation strength? (b) Is this 2D transition con¬ 
tinuous and critical (as in ID), or first order? (c) Does 
the spontaneous misalignment affect the transition, and 
when so, to what effect? 

Here we present simulations for a model incommen¬ 
surate monolayer which show that: (a) Super lubricity is 
replaced by pinning through a sharp transition, (b) That 
transition is, at least for the set of parameters used, of 
first order rather than continuous, (c) Spontaneous rota¬ 
tions very definitely affect the transition point, decreas¬ 
ing its critical value of corrugation Uq. 

The paper is organized as follows. In Sect. [H] we intro¬ 
duce the model and the protocols adopted in the molecu¬ 
lar dynamics simulations. In Sect. Hill we present results 
which indicate a first-order structural transition taking 
place in the colloid monolayer as a function of increas¬ 
ing corrugation. Section |TV] discusses the coincidence be¬ 
tween this structural transition and the first-order onset 
of a static friction force Fg. In Sect. El we describe the 
dependence of the static properties of the monolayer on 
the misalignment angle. Section IVTl contains a discussion 
of possible experimental verifications, as well as our final 
remarks and conclusions. 


II. MODEL AND SIMULATIONS 

We describe the colloidal dynamics using the same 
model and parameters as in Refs. 0 and M In short, the 
2D layer of colloidal particles is represented by screened- 
Coulomb repulsive, classical point-like particles, moving 
in 2D in a periodic potential (the corrugation) of ampli¬ 
tude Uq, periodicity oi, and triangular symmetry. The 
overdamped dynamics of these particles is generated by 
integrating T = 0 damped equations of motion, with a 
large viscous damping coefficient rj = 28. All results 
are expressed in terms of the system of units defined in 
Table I of Ref. M Simulations were performed adopt¬ 
ing 2D periodic boundary conditions (PBC). Any de¬ 
sired misalignment angle between the corrugation and 
the colloidal lattice is implemented by means of a suit¬ 
ably chosen supercell, built using a standard procedure 
as follows.— The two lattices are defined by the pairs 
of primitive vectors ai = ai(I,0), &2 = ai(0.5,-x/S), and 
bi = ac(cos0, sin6>), b 2 = ac(cos(0 -|- tt/S), sin(0 -|- tt/S)). 
The overall supercell-periodic structure can be realized if 
four integers exist which satisfy the matching condition 
niai + n 2 SL 2 = mibi -|- m 2 b 2 . The ensuing supercell is 


(ni,n2) 

{mi,m 2 ) 

P 

e 

Np 

(96,96) 

(89,89) 

0.92708333 

0 

7,921 

(105,59) 

(107,71) 

0.92705245 

2.54° 

20,701 

(96,72) 

(100,54) 

0.92705012 

5.07° 

18,316 

(78,74) 

(59,103) 

0.92705113 

9.78° 

17,332 


TABLE 1. Parameters of the four supercells adopted for the 
simulations in PBC (see text for definitions). The supercells 
are determined following the procedure outlined in the text. 
For 9 = 0, the adopted mismatch p = axja^ = 89/96 corre¬ 
sponds to the sixth approximation in the continued fraction 
expansion of 3/(1 -I- v^). 


that of a triangular lattice of size L = Imibi -|- m 2 b 2 |, 
containing a total number of particles Ap = m 2 -l-mim 2 -l- 
m\. In practice we fix Oc = 1 and vary ni_ 2 , 'mi ,2 in search 
of structures with a mismatch p = ai/^~3/(I-|- v^) - 
close to the experimental values of Ref. Ij - and 9 near the 
desired value, with the obvious additional constraint that 
the number of particles Np should not be too large. We 
consider in practice the aligned configuration 0 = 0, plus 
the following misaligned configurations: 0 ~ 5°, 6* ~ 10° 
and 9opt ~ 2.54°. The latter is close to the (Novaco- 
McTague) equilibrium misalignment angle 0 nm — 2.58° 
predicted by weak-coupling elastic theory— Table |T] 
collects the adopted supercell parameters. 


As mentioned earlier, the mismatch in this case pro¬ 
duces a moire pattern which could be roughly described 
as a patchwork or superlattice of small domains where 
colloids and corrugation are mutually (nearly) commen¬ 
surate, separated by a network of anti-soliton lines whose 
thickness decreases with increasing corrugation strength. 
From the experimental point of view, the underdense 
regime chosen here for exemplification is better suited 
than the overdense case p > 1 where local compressions 
may lead to buckling of particles out of the plane, and 
where the energetics (not symmetrical with respect to 
p > I) is less convenient. The initial step of our study 
is to identify the distorted colloid lattice configuration of 
lowest energy in presence of the corrugation. For that, 
we carry out T = 0 damped-dynamics simulations test¬ 
ing different protocols. In the first we increase the cor¬ 
rugation amplitude Uq in steps AU = 0.018 - eventually 
reduced near transition critical points - then reversing 
the sign of AU and decreasing Uq back to zero in or¬ 
der to check for hysteresis across structural transitions. 
The maximum values of the corrugation considered is 
ffmax = 1. In the second protocol we switch on abruptly 
the corrugation at given values of Uq with the monolayer 
initially either undistorted or in a configuration previ¬ 
ously fully relaxed at I/max- For each Uq we eventually 
take the lowest-energy configuration. Eventually the best 
configurations were always obtained using the first pro¬ 
tocol described above. 














FIG. 2. (Color online) Examples of the moire patterns ob¬ 
tained at p « 0.927 for misfit angles (a) 6 = 0, (b) 0opt ~ 
2.54°, (c) 0 ~ 5°, and (d) 6 ~ 10°. A small portion of the 
simnlation snpercells is displayed, containing an undistorted 
monolayer. The particles’ colors reflect the underlying corru¬ 
gation potential W{x,y)'. dark for potential minima, bright 
for maxima. According to Eq. m, at 0 ~ 2.54° the moire 
orientation is V’ — 30°. As 6 increases the superstructnre pe¬ 
riodicity shrinks rapidly and rotates all the way to V’ — 60° 
(c). 


The depinning of the monolayer is studied by applying 
to each colloid a driving force Fd, generally along a high- 
symmetry direction of the laser substrate potential. For a 
given value oi Uq, the single particle barrier for the onset 
of motion is Fgi = SirUo/Oai. The external force is then 
increased adiabatically in steps AF = 0.002, much lower 
than the lowest value of Fsi considered. The duration is 
fixed in such a way that a single free particle would slide 
by a total distance of 10 Oc. A colloid configuration is 
classified as sliding, when the procedure produces a final 
center-of-mass displacement Aa;cm ^ 2ac. 

Preliminary to presenting our results, it is useful to 
recall a few geometrical concepts. A given mismatch p 
and misalignment angle 6 induce a specific moire pattern 
of the colloidal crystal. The moire superstructure of an 
undistorted monolayer is entirely described by its overall 
periodicity L and its orientation angle '0 relative to the 
substrate. Similar to the beats of two sinusoids in one 
dimension, the values of 0 and L are governed by the 
difference G — G' between the smallest reciprocal lattice 
vectors of the monolayer and optical potential. The rela¬ 
tions connecting p and 6 to L and 0 are given for example 
in Ref. [U Here we recall the angular relation 


cos 6 = p ^ sin^ 0 -I- cos 00 1 — p~‘^ sin^ 0 , (1) 

which, for p ~ 1 produces a rapid reorientation of the 
moire pattern for small changes of 0 ~ 0, as visualized in 



FIG. 3. (Golor online) The total energy (Etot) per particle of 
the colloidal monolayer optimized in PEG, for a few misalign¬ 
ment angles 6. Gurves correspond to increasing corrugation 
amplitude Uo ~ 0.1 — 0.3. Energies are measured relative to 
that of the perfect crystal at Uo = 0, and are scaled with 
respect to Uo for better comparison. 


Fig. [21 


III. RESULTS: THE TRANSITION IS 
STRUCTURAL 

We begin by describing the evolution of the ground 
state configuration of our model underdense colloidal 
monolayer as a function of increasing corrugation 
strengths Uq, obtained at fixed incommensurability p ~ 
0.927 and for a grid of angular misalignment angles 9. 
Figure [3] shows the total potential energy per particle 
(Etot) as a function of Uq and for several small 6 val¬ 
ues. The lowest energy is seen to occur at a small but 
finite 0opt- This is a small but significant misalignment 
angle fairly close to the theoretical, weak-coupling mis¬ 
alignment 0 nm suggested by the theory of orientational 
epitaxyjiSi^— The detailed angular dependence of total 
colloid energy in the optical lattice was investigated in 
Ref. [H. The mainly longitudinal deformation for 6 = 0, 
is replaced in the optimal misalignment 0opt, by ener¬ 
getically cheaper largely transverse deformations. This 
optimal misalignment angle is 0opt = 2.54°, to a good 
approximation independent of corrugation Uq. 

We wish to understand now what happens to the col¬ 
loid structure, energy, and eventually to the static fric¬ 
tion as a function of Uq at 9 = 0opt. For that purpose, 
we consider that the total potential energy per particle 
(Etot) = Etot/Np is composed of two separate contribu¬ 
tions: 

{Etot) = {W) + (Upp), (2) 

where (FF) is the average colloid-optical lattice interac¬ 
tion energy, and {Upp) the average colloid-colloid repul¬ 
sion. Figured^ displays the two contributions to the to¬ 
tal energy {Etot)- The total energy evolution is smooth 
until Uq = 0.22 where we observe a sudden decrease in 
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FIG. 4. (Color Online) (a) The total potential energy (Etot) 
per particle, its periodic-potential (corrugation) contribution 
{W), and interparticle contribution (I/pp) as a function of the 
corrugation Uq for the optimally misaligned colloidal mono- 
layer {6 = 6opt = 2.54°). The sudden drop of {W) and the cor¬ 
responding jump of {Upp) signal a first-order transition, tak¬ 
ing place at constant incommensurability and misalignment 
angle. As will be shown later, the weak-corrugation phase 
is unpinned and superlubric, the strong-corrugation phase is 
pinned by a finite static friction, (b) The intersection be¬ 
tween the total energy branches representing the superlubric 
(small Uq) and statically pinned (large Uq) structures. Like in 
all first-order transitions, each phase survives in metastable 
state beyond the transition point. 


{W) accompanied by an increase of the repulsive term. 
This indicates a first-order structural transition in the 
monolayer. Figure [S^ displays the 0opt moire patterns 
of the relaxed monolayer for Uq = 0 . 108 , a weak cor¬ 
rugation well below the first-order transition. For this 
weak corrugation the triangular colloidal crystal is only 
slightly affected by the substrate potential. Panel [SJd il¬ 
lustrates the corresponding colloid displacements relative 
to the rigid unperturbed crystal (Uq = 0 ). The misalign¬ 
ment of the colloidal crystal over the lattice corrugation 
is only slightly modulated by a vortex-like distortion re¬ 
flecting the shear nature of the solitons, whose width is 
larger than their separation. Panel[SJ: portrays the moire 
pattern for Uq = 0 . 395 , a strong corrugation well above 
the transition. A hexagonal superlattice of domains of 
colloids locally commensurate with the corrugation are 
separated by “atomically” sharp antisolitons lines. 

It is also instructive to look at the average local orien- 



FIG. 5. (Golor online) The moire patterns of the fully-relaxed 
incommensurate colloid monolayer at 0opt = 2.54°. A small 
portion of the simulation supercell is shown. Each dot rep¬ 
resents a particle, colored according to the local corruga¬ 
tion potential W{x,y): dark for potential minima, bright 
for maxima, (a) The optimal weak-corrugation configuration 
{Uo = 0.108). (b) Zoomed-in region with arrows highlighting 
the particles’ displacements (magnified 20 x) from the ideal 
triangular lattice to the fully-relaxed configuration. Vortex¬ 
like displacements tend to locally rotate the domains back 
into registry with the substrate, (c) The optimal strong- 
corrugation state {Uq — 0.395). The hexagonal domains ap¬ 
pear now fully commensurate, and separated by sharp domain 
walls. Here most colloids fall very close to a potential mini¬ 
mum. 
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FIG. 6. (Golor online) The average local orientation of the 
monolayer Oioc (see text for definition) plotted as a function 
of corrugation. The vertical dashed line separates the weak- 
corrugation and the strong-corrugation phases. 


tation of the monolayer (see Fig. [6]), defined as 

(hi) 

where the sum is over all M pairs {i,j) of nearest neigh¬ 
bor particles, and 9ij is the angle between the relative 
position vector and the j;-axis. At the transition, 

the commensurate domains realign significantly with the 
triangular potential. Figure IHd shows a crossing of the 
total energies of the weak-corrugation and the strong- 
corrugation phase at the transition point, as a function of 
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FIG. 7. (Color online) The static friction force Fs of the p ~ 
0.927 colloidal monolayer, normalized to the single-particle 
barrier Tsi. The dashed line indicates the transition point 
corrugation Uc = 0.22 where the colloid monolayer suddenly 
locks onto the corrugation. The error bar is defined by the 
force increment AF = 0.002 used in the simulation protocol. 


Uq. a separate stability analysis shows that each phase 
survives as a metastable state over a finite corrugation 
interval on the “wrong” side of Uc, eventually collapsing 
at its respective spinodal point. 


IV. RESULTS: SUPERLUBRICITY AND 
STATIC FRICTION 

We now come to the static friction force of the col¬ 
loid monolayer, the main interest of the present investi¬ 
gation. It is obvious that the first-order structural tran¬ 
sition uncovered in the previous section must affect it. 
Figure 0 shows the static friction force Fg obtained by 
simulation with p ~ 0.927 and 0opt = 2.54° as a function 
of the corrugation strength Uq. For weak corrugation 
Uq < Uc the static friction is small below the error. In 
this regime the colloid monolayer slides freely already at 
the smallest applied force, demonstrating that the incom¬ 
mensurate state is indeed superlubric^ At the transition 
corrugation the static friction suddenly jumps to O.OSFgi. 
The monolayer becomes pinned, and the superlubricity 
is lost through a first-order transition. Unlike ID, where 
the pinning transition was of second order, here it coin¬ 
cides with the first-order structural transition. 

Leaving the transition order aside, it is easy to show 
that the underlying reason for pinning in the present 2D 
case is still quite similar to that in the ID FK model, 
namely a collapse of the dynamically-available phase 
space. In ID, that collapse was signaled by the simultane¬ 
ous collapse of the probability for any particle to occupy 
a local potential maximum. In the present 2D case we 
can similarly investigate the probability for colloids to oc¬ 
cupy, in their optimal configuration, the {x, y) triangular 
region above the saddle-point among adiacent local min¬ 
ima of the corrugation potential. If the analogy with ID 
holds, a superlubric, ergodic state should populate abun¬ 



FIG. 8. (Golor online) (a) The fraction of colloidal parti¬ 
cles which sample repulsive regions of the substrate potential, 
i.e. those positions where W(x,y) exceeds the saddle energy, 
plotted as a function of the corrugation strength Uq. (b) 
The triangular potential W{x,y). Dark and bright regions 
correspond respectively to minima and maxima, isolines are 
reported at the saddle-point value and a primitive cell is high¬ 
lighted at the center of the density plot, (c-f) The positions of 
the particles reported inside one primitive cell of the substrate 
potential. Gonfigurations are shown for different corrugation 
values: (c) Uq = 0.108 - superlubric phase; (d) Uq = 0.219 
- just below the transition; (e) Uq = 0.221 - just above the 
transition; (f) Uq = 0.395 - well above the transition. 


dantly this region, while in a pinned, broken-ergodicity 
state, the population probabilty at the center of these 
regions should collapse to zero, and even remain zero 
across a whole neighborhood of the maximum - a finite 
’’disorder parameter” Indeed, as shown in Fig. [5^, at 
the first-order superlubric-pinned transition, these repul¬ 
sive regions around the maxima depopulate dramatically. 
This is represented pictorially in Fig. [Bt-f, reporting the 
particles’ positions folded in one primitive cell of the cor¬ 
rugation lattice (see Fig.|8}3). These properties represent 
a 2D version of the disorder parameter If, which mea¬ 
sures the width of the largest gap appearing in the Hull 
function at the Aubry transition in the incommensurate 
ID FK model4 It should be noted that in the pinned 
colloid state the static friction force remains well below 
the single-colloid limit value Fgi, as shown in Fig. [T] The 
reason for this lower static friction is that of course the 
monolayer sliding is still strongly collective, whereby the 
particles surrounding a given particle attempting to cross 
a forbidden region actually help by pushing it across the 
barrier. 
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FIG. 9. (Color online) Unrotated 0 = 0 mismatched col¬ 
loid monolayer, p « 0.927. (a) Fraction of colloidal particles 
which populate repulsive regions of the corrugation potential 
(where W{x,y) exceeds the saddle-point energy) plotted as a 
function of the corrugation strength t/o across the superlubric- 
pinned transition. The dashed line marks the transition point 
Uc = 0.391. (b) Static friction force Fs, normalized to the 
depinning force Fsi of a single particle in the same periodic 
potential. Above the transition the monolayer locks onto the 
corrugation. The error bar is dehned by the force increment 
AF = 0.001 used in the simulation protocol. 


V. RESULTS: STATIC FRICTION AT 
DIFFERENT MISALIGNMENTS 

An interesting aspect of the transition to the 2D locked 
state is its dependence on the misfit angle 9. As discussed 
in Sect. uni at the optimal orientation 0opt the colloidal 
crystal deforms mainly via the softer transverse modes, 
facilitating the interdigitation with the substrate. If the 
monolayer is rotated away from 0opt its grip of the corru¬ 
gation weakens: as a result the superlubric state survives 
up to a larger corrugation strength C/q. As the simplest 
example, consider the aligned monolayer configurations 
obtained at 0 = 0. 

Figure shows the evolution of the number of parti¬ 
cles sampling repulsive regions of the corrugation poten¬ 
tial. In this case too we identify a corrugation Uc = 0.391 
above which the monolayer pins, again with a first-order 
transition. As expected, the transition corrugation is 
larger - almost a factor of two - than that at the op¬ 
timal misalignment, Uc = 0.22. Figure [SJd shows the 
static friction force Fg across the structural transition. 
The aligned monolayer is superlubric up to Uc, where we 
observe a sudden surge of the static friction Fg . Actually 
in this metastable aligned case the subsequent evolution 
of static friction becomes somewhat erratic. This be¬ 
haviour can be explained as follows. A proliferation of 
energetically close metastable configurations - each char¬ 
acterized by its own static friction ~ is expected^ above 
the transition point. While at 0opt all protocols produced 
configurations with regularly-spaced hexagonal domains. 



FIG. 10. (Color online) Panels (a-d): the moire patterns of 
the fully relaxed monolayer at 0 = 0. A small portion of the 
simulation supercell is shown. Each dot represents a particle, 
colored according to the local corrugation potential W{x,y)\ 
dark for potential minima, bright for maxima. Panels (e-h): 
the particles’ positions are reported into one primitive cell of 
the substrate triangular potential. Gonhgurations are shown 
for the following values of the corrugation: (a,e) f/o = 0.108 
- superlubric phase. (b,f) Uo = 0.386 - just below the transi¬ 
tion. (c,g) Uo = 0.395 - just above the transition, where the 
energetically favored configuration consists of an arrangement 
of irregular domains, (d-h) Uo = 0.712 - a strong-corrugation 
phase corresponding to an arrangement of regularly-shaped 
registered domains. 


at 0 = 0 we observe several both irregular and regular 
arrangements of registered domains, which then happen 
to be energetically favored at different corrugation values 
Uq (see Fig.fTUh.d). Some fluctuations in Fg are therefore 
expected when one also considers that in incommensurate 
systems with many degrees of freedom, states with lower 
energy do not necessarily display a higher static friction 
forcei^ This complexity disappears below Uc where there 
always is a unique structure, shown in Fig.fTOh.b. respec¬ 
tively away and close to the transition. Summing up this 
part, the evolution of the monolayer structure as a func¬ 
tion of C/q is highly sensitive to the misalignment angle. 






















VI. DISCUSSION AND CONCLUSIONS 


Colloids in optical lattices are systems recently in¬ 
troduced, potentially representing powerful tools which 
may contribute to our understanding of friction between 
crystalline surfacesColloidal monolayers in an in¬ 
commensurate lattice corrugation potential are shown to 
represent ideal systems to study the superlubric-pinned 
transition, expected as a function of increasing corruga¬ 
tion strength, but never realized so far in any 2D sys¬ 
tem, under constant geometrical conditions. In ID it is 
well known that at this transition the onset of a disor¬ 
der parameter - measuring a region of phase space inac¬ 
cessible to particle positions - impedes the dynamics of 
free-sliding motion, giving rise to static friction through a 
continuous phase transitional^ In the 2D case of colloids 
we predict a transition of similar character, the main dif¬ 
ference being now its first-order character. The relative 
alignment and misalignment of the two lattices plays an 
additional role. We found vanishingly small static friction 
up to a critical value Uc corresponding to a joint struc¬ 
tural and dynamical transition of the monolayer. At the 
transition the colloidal crystal adapts to the periodic po¬ 
tential by a local rotation (back to registry) of the nearly 
commensurate domains forming the moire of the incom¬ 
mensurate phase. 

Similarly to the ID case, pinning of the monolayer is 
due to a collapse to zero of the number of particles sam¬ 


pling the most repulsive regions of the periodic poten¬ 
tial. These poorly placed particles are those which are 
easily set into motion in the superlubric phase. Their 
disappearance is therefore responsible for the arisal of 
static friction in the pinned state. Experimentally, van¬ 
ishingly small values of Fg have been observed in 2D 
incommensurate structures,— but the realization and 
the analysis of the transition between superlubric and 
pinned states is missing to our knowledge. The present 
work suggests the possibility to pursue and investigate 
in detail the superlubric-pinned transition by gradually 
increasing the optical-lattice strength. Additional phe¬ 
nomena may be accessed by forcing rotations and mod¬ 
ifying the natural misalignment of the colloidal crystal 
with respect to the optical lattice. Since the superlubric 
phase at 0 7 ^ 0opt survives up to larger corrugation val¬ 
ues, rotations of mismatched sliders initially away from 
the optimal orientation might in fact end up in reentrant 
pinning as soon as 0opt is hit, in a way similar to the 
matched case (p = 1, 0opt = 0) of a rotated graphene 
flake sliding over graphite<^ 
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